Mergers of Magnetized Neutron Stars with Spinning Black Holes: Disruption, 

Accretion and Fallback 
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We study the merger of a neutron star in orbit about a spinning black hole in general relativity 
with a mass ratio of 5 : 1, allowing the star to have an initial magnetization of 10^^ Gauss. We present 
the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. We 
see no significant dynamical effects in the simulations or in the gravitational waveform resulting 
from the magnetization. We find that only a negligible amount of matter becomes unbound; 99% 
of the neutron star material has a fallback time below 10 seconds to reach the region of the central 
engine and that 99.99% of the star interacts with the central disk and black hole within 3 hours. 



Introduction. The spectacular energetics associated 
with short gamma ray bursts (sGRBs) are difficult to 
explain, requiring complex models synthesizing a vari- 
ety of different components (see e.g. [1^3]). Key among 
these is the inclusion of extreme gravity responsible for 
accelerating plasma to high Lorentz factors. Gonsensus 
is building for a scenario in which the gravitational field 
results from the merger of two highly compact objects: 
either a black hole and a neutron star (BH-NS) or a bi- 
nary neutron star system (NS-NS). These systems radiate 
strongly in electromagnetic and gravitational wave bands 
making them ideal candidates for combined observations 
(e.g. y). The validation of such models requires a care- 
ful comparison of both electromagnetic and gravitational 
wave signatures with theoretical predictions. 

Such predictions require sophisticated simulations in- 
corporating the necessary physics ingredients. A the min- 
imum, they require solving the full, nonlinear field equa- 
tions of general relativity along with relativistic hydro- 
dynamics. For the particular case of BH-NS, numerical 
models have recently begun achieving interesting success, 
and, despite the complexity of the parameter space in- 
volved, a common picture is emerging towards connect- 
ing the system with sGRBs. For instance, recent results 
indicate that the disk resulting from the merger of a 
non-spinning BH with a NS (approximated by a F-law 
ideal fluid) is far less massive than what leading sGRB 
models require On the other hand, if the BH is 

(sufficiently highly) spinning, the resulting disk is signif- 
icantly more massive and falls within values consistent 
with the leading sGRB models resulting from BH-NS 
merger (e.g. [lol-llij: see [l3|,|l3 for recent reviews on the 
subject). Gomplementary efforts to understand possible 
observable electromagnetic counterparts are actively be- 
ing investigated (e.g. |15l4l7j). 



Beyond the importance for their connection to sGRBs, 
BH-NSs are also one of the most likely sources of de- 
tectable gravitational waves with current and near-future 
earth-based gravitational wave detectors. Similar to bi- 
nary black hole coalescence, BH-NS mergers (with a 
black hole mass above 10 — 20 Mq) are very bright 
gravitational wave sources, but they are also expected to 
demonstrate remarkable sensitivity to the details of the 
neutron star due to tidal effects within the most sensitive 
frequency window of these detectors (e.g. [ill, fisl. [l9| ) . 

These systems are complex with a diverse phenomenol- 
ogy. Indeed, the equation of state of the fluid, the 
spin of the BH, a nonvanishing magnetic field, and 
neutrino cooling all can have a profound influence on 
the dynamics of the system. Different studies have 
been presented to explore the phenomenology related 
to the first two points above. In the current work 
we further explore these options and study the possi- 
ble impact of the star's magnetic field throughout the 
merger with a BH, employing our General Relativistic- 
MagnetoHyDrodynamics code (GRMHD). This allow us 
to study the system with effects that dominate the dy- 
namics (general relativity and hydrodynamics) together 
with magnetic effects which may play a role through the 
merger and early-post merger. We explore the dynamics 
of material disrupted from the star. In particular, we 
estimate typical fall back times, and consider whether 
the observed dynamics is consistent with processes sug- 
gested as drivers of sustained emissions after the main 
burst from sGRBs (e.g. [13, lloHill ) . 

We model the neutron star material using relativis- 
tic ideal MHD, coupled to the full Einstein equations 
of general relativity to accurately represent the strong 
gravitational effects during the merger. Our numerical 
techniques for solving these coupled equations have been 
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(a) t=5.9 ms (b) t=16.3 ms (c) t=27.3 ms 




FIG. 1: Isosurfaces of density (6.18 x 10^*^ g/cm^) and appar- 
ent horizon (black spheroid) at various times for the magne- 
tized evolution. 



thoroughly described and tested previously 23-27|. 
Set-up. We consider a BH-NS system where the neutron 
star is possibly magnetized and described by F = 2. We 
begin with quasi-circular initial data constructed with 
LORENE [28]. We adopt a realistic mass ratio [29| q = 
Mns/^bh = 1/5, and, for the magnetized cases, we add 
an initial, poloidal magnetic field to the neutron star by 
assuming a purely azimuthal vector potential as A^^ = 
tz7^ max(P — Pvac, 0), [30| with w the cylindrical radius 
and pressure Pyac/c^ 10^ g/cm^. This yields a field 
confined to the stellar interior with maximum magnitude 
of 10^2 G. The BH has spin 0.5 and mass Mbh = 7Mq. 
The neutron star baryon (gravitational) mass is 1.473 
(1.334) Mq). The total mass of the system (BH mass 
plus NS gravitational mass) is Mt = Mbh + Mns = 
8.33M0 and compactness (M/R) is 0.1. 

The initial data are evolved in a cubical computational 
domain defined by G [—443 km, 443 km], and we em- 
ploy adaptive mesh refinement which tracks the two com- 
pact objects. However, once the star disrupts, the fluid 
is no longer particularly localized which constrains the 
size of the grid spacing of the coarsest level. We adopt 
a coarse grid with spacing of A = 2.952km which cov- 
ers the entire computational domain. Higher resolution 
is achieved by adopting 3 further levels of refinement for 
which the finest spacing has A = 0.738 km (convergence 
comparisons were made with respect to runs with just 1 
and 2 levels of refinement). The refinement criterion is a 
combination of the gradient of the metric, the amplitude 
of extrinsic curvature, and the density of the fluid. 
Results. We focus primarily on the case in which the 
BH has spin aligned with the orbital angular momentum. 
We note that the main features discussed next are essen- 
tially the same for both magnetized and unmagnetized 
cases except at low density values; thus, unless noted the 
results discussed above stand for both cases. 

For our spinning BH configuration, estimates of the lo- 
cations of the inner-most stable circular orbit (ISGO) [HI, 
Isij and the mass-shedding limit [111 indicate that 
both are at comparable distances from the BH (:^ 
(3.36, 2.86)Mt respectively). Thus, the neutron star is 
expected to disrupt around the time it crosses the ISCO. 
Were the BH not spinning, as discussed previously, the 
NS generally falls into the BH quite quickly with little 



mass remaining in any accretion disk. 

As shown in Fig. [H the neutron star orbits the BH for 
about two orbits and tidal disruption begins at time t 9 
milliseconds. We extract the resulting gravitational wave 
signal by computing the Newman-Penrose Weyl scalar 
ip4 at different coordinate radii and then decomposing 
onto an appropriate spin- weight —2 basis. The dominant 
mode of this signal is illustrated in Fig. [2fleft). Until 
about :^ 12 ms when the star approaches the ISCO, one 
sees a signal generally characteristic of the quadrupole 
radiation of two orbiting masses. However, the star then 
crosses the ISCO and soon after begins to shed. This 
leads to a rapid decrease in the gravitational wave out- 
put [5|,[lo|,[33|. Notice that the familiar ringdown pattern 
observed in binary black hole mergers is essentially ab- 
sent due to the continuous in-fall of material. 

Figure [Sfright) displays the power spectrum of the 
gravitational wave strain. Also shown with vertical bars 
are estimates of two frequencies, /[^qq and /qnm which 
characterize the system. These are obtained via simple 
first-principles estimates based on orbital frequencies cor- 
responding to ISCO and "light-ring" locations or by ex- 
amining the obtained solution. As discussed in [l9|, [sij a 
simple estimate can be obtained by an "angular momen- 
tum balance" argument at the ISCO for the two-body 
problem, ignoring radiative and disruption effects. This 
estimate provides a value for the final BH spin of ~ 0.7 
which gives /[^qq — llOOHz. An accurate number, on 
the other hand, can be obtained by a direct inspection of 
the horizon at t=27.3 ms which indicates a ratio of po- 
lar to equatorial circumference of the BH of 0.931. This 
ratio corresponds to a BH with spin a/M = 0.56 which 
would indicate /jg^o — 880Hz. We thus expect a quali- 
tative (smooth) change in the wave strain in this range 
of frequencies as the system transitions from an orbit- 
ing pair to an accreting BH, and indeed our spectrum 
shows such a change as the overall slopes before and af- 
ter these frequencies are markedly different (~ —1/6 vs 
~ —3.5). (Notice however early oscillations are evident 
due to both eccentricity from the initial configuration 
and neutron star oscillations) Similarly, the estimates of 
/qnm — for the dominant quasi-normal frequency — for a 
BH with spin a/M = (0.56,0.7) are (1900, 2100)Hz. 

Additionally, as the system is asymmetric, there is a 
net flux of momentum carried out by the gravitational 
waves, which results in the final hole acquiring a recoil 
velocity. A 2PN estimate of the kick [s^ to 900Hz gives 

2km/s which is expected as this value does not take 
into account the merger stage. The actual value com- 
puted from the extracted waveform is 40 km/s which 
is below 150km/s that would be estimated for a binary 
black hole system with otherwise equal physical parame- 
ters [35|. That the BH-NS recoil value obtained is lower 
than the analogous BH-BH system is expected as the for- 
mer radiates less energy and momentum than the latter. 

A significant amount of the NS matter is accreted 
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FIG. 2: (left) The / = 2,m = 2 mode of r^4 for mergers 
examined here (removing the initial artificial stage related to 
the initial data). Both the magnetic and nonmagnetic wave- 
forms are the same; extraction performed at coordinate dis- 
tances 295 km, 332 km, and 369 km and adjusted for travel 
time, (right) Power spectrum of the wave-strain f\h\ in which 
in grey bands denote frequencies associated with the ISCO 
and quasi-normal ringing corresponding to a BH with mass 
Mt = 8.33M0 and spin of a/Mx = (0.56, 0.7). 



through the merger process, and the spin of the BH 
consequently increases. As mentioned, the horizon ge- 
ometry indicates a BH with spin a/Mx = 0.56. This 
value is consistent with the analogous case for a binary 
BH merger for which a simple model predicts a final BH 
spin of a = 0.7 M [3l|. That the value here is smaller 
is expected because the angular momentum of the fluid 
remaining outside the BH is not captured. 

As a result of the merger, for the spinning cases a sig- 
nificant amount of matter (O.ITM©) at about t 20ms 
is observed, regardless of the magnetization considered. 
Of this, we can estimate that a disk is formed with a 
mass equal to about 1% of the initial stellar mass for 
both the magnetic and non-magnetic cases (based on 
the integrated fluid mass on the finest resolution mesh 
about the BH). However, we note that significantly more 
mass than this, about O.OTM©, remains outside the BH 
as shown for the magnetized case in Fig. [3] for late times 
(times t 30ms). The vast majority of this material 
has remained gravitationally bound to the BH forming a 
reservoir that will eventually return to interact with the 
central engine (see also [l2|). 

It is useful to examine more closely the disk structure, 
temperature and velocity profile. We find that the struc- 
ture consists of a hot and vertically thick region where 
material in the spiral arm has shocked due to intersection 
of stream lines while much of the tidal debris remains thin 
and cold prior to shocking. The temperature, estimated 
from the ratio of pressure to density, varies between 10^^ 
- 10^^ K while the tidal tail of material thrown off is 
substantially cooler, around ~ 10^ K. The velocity pro- 
file of the disk and tidal tail in both magnetic and non- 
magnetic cases is shown in Fig. |4j As the merger and 
subsequent disk formation proceeds the magnetic field 
is redistributed into a toroidal configuration and grows 
through magnetic winding. Figure 4a illustrates the mag- 
netic field structure at t = 22.2ms. Furthermore, as a 
result of the disruption and merger process, a significant 



amount of matter is thrown out of the region close to the 
BH but most all of it (greater than 99%) remains bound 
as its speed is below the escape speed of the BH. This 
bound material will eventually interact with the central 
engine again. We calculate a fallback time for individ- 
ual fluid elements based on the method detailed in (ioj . 
Fig. [5] shows distributions of the disrupted matter for a 
few times during the evolution (tidal disruption of the NS 
has begun by 8.9 ms while the spiral arm of tidal debris 
is well- formed at 16.3 ms). These distributions indicate 
that the accreted mass follows a power-law in fallback 
time such that the accretion rate falls off with exponent 
about —5/3; we find a similar exponent in the case of a 
non- spinning BH as well. This value equals that of Phin- 
ney [36| for accretion of material stripped from a main 
sequence star by a supermassive BH This agreement is 
quite interesting due to a number of differences between 
the system studied in [36| and here: (i) our estimates de- 
rive from a relativistic evolution; (ii) our BH is spinning; 

(iii) our system is in a quasi-circular orbit, not parabolic; 

(iv) our disrupted star is a NS, not a main sequence star 
and (v) the mass-ratio is far closer to unity such that the 
physical sizes of the BH and NS are comparable. 
Final comments. Our results^together with other 
studies of BH-NS mergers [1, [IqI, , indicate that BH- 



NS systems with realistic mass-ratios give rise to a suf- 
ficiently massive disk for connecting with sGRBs if the 
spin of the BH is sufficiently high (otherwise the mass 
in the resulting disk decreases considerably). For these 
cases, gravitational waves from the system will manifest 
subtle differences in the waveforms tied to the equation 
of state of the star as the mass-shedding radius is not 
far from the ISCO. Detecting such differences, however, 
will require delicate work on the data analysis front as 
the frequency window in which these differences will arise 
is quite small. This issue is also encountered in binary 
neutron star systems [33|. Furthermore, the observation 
that the burst in gravitational waves might not be fol- 
lowed by quasinormal black hole ringing bears relevance 
to burst searches in data analysis pipelines adopting dif- 

33, la 



ferent models to capture the burst behavior 

It is interesting to consider these results in the con- 
text of short GRBS. Our evolutions of BH-NS merg- 
ers indicate several interesting stages: (i) at early times 
in the merger {c^ 10~^ — 10~^s), the BH hyperaccretes 
suggesting it might be a good candidate for creating 
a fireball through neutrino annihilation based on the 
mass accretion rate and remnant mass [39|; (ii) at later 
times (< lO^s), sufficient mass will be falling back which 
might support long sustained emissions via r-processes, 
consistent with observed emissions in roughly 30% of 
short GRBs ^|; (iii) at even later times (> lO^s up to 
about lO^s), there remains enough bound mass (roughly 
10~'^Mq) to be consistent with estimates of electromag- 
netic merger counterparts to gravitational waves [l7j . 
Acknowledgments. We dedicate this work to the 
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FIG. 3: The integrated mass (solid line) vs. time for the 
magnetic case. The plus signs indicate the integral of material 
only outside the ISCO. The dashed curve shows the integral 
only over material moving slower than the escape speed for a 
7 solar mass BH. The dot-dashed line (lowest curve) indicates 
unbound material (moving at or faster than the escape speed, 
regardless of direction) multiplied by a factor of 4 to put it on 
the same scale. The equivalent integrals for the unmagnetized 
simulation yield nearly identical values up to t :^ 20ms. 
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FIG. 4: System behavior at t=22.2 ms (axis in kms and den- 
sity color key in units of g/cm^). (a) density isosurface and 
magnetic field lines with the BH indicated by the central black 
spheroid, (b) - (d) fluid density grouped according to the 
speed of the fluid in the equatorial plane. There is no fluid 
with velocity more than 0.5. 
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FIG. 5: Fallback accretion histogram. For all bound mass 
beyond 1.5 times the ISCO, a fallback time is estimated and 
the results are binned as shown. Distributions are computed 
for three times shown by the different colors that roughly 
correspond to the start of the tidal disruption event, the peak 
of gravitational radiation emitted from the system and the 
point where a disk encircles the central BH. The dashed line 
indicates a power-law with a slope of —2/3 for reference. 
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